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Abstract 

Using a multidimensional cut-off technique, we obtain expressions for the cut-off dependent 
part of the vacuum energy for parallelepiped geometries in any spatial dimension d. The 
cut-off part yields nonrenormalizable hypersurface divergences and we show explicitly that 
they cancel in the Casimir piston scenario in all dimensions. We obtain two different 
expressions for the d-dimensional Casimir force on the piston where one expression is more 
convenient to use when the plate separation a is large and the other when a is small (a useful 
a — ► 1/a duality). The Casimir force on the piston is found to be attractive (negative) for 
any dimension d. We apply the (i-dimensional formulas (both expressions) to the two and 
three-dimensional Casimir piston with Neumann boundary conditions. The 3D Neumann 
results are in numerical agreement with those recently derived in arXiv:0705.0139 using 
an optical path technique providing an independent confirmation of our multidimensional 
approach. We limit our study to massless scalar fields. 
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1 Introduction 



With idealized boundary conditions, such as Dirichlet and Neumann boundary conditions, vac- 
uum energy calculations have cut-off dependent terms that diverge in the limit as the cut-off 
A tends to infinity. These divergences can be classified as either bulk (volume) divergences 
or lower-dimensional surface divergences (for simplicity, all divergences besides the volume 
divergence will be classified as hypersurface or simply surface divergences regardless of the di- 
mension). The volume divergence poses no problems since the Casimir energy, by definition, is 
obtained after subtracting out the volume divergence from the vacuum energy. In other words, 
the volume divergence is renormalizable and can be set to zero by simply adding a constant 
counterterm to the Hamiltonian. In contrast, the surface divergences are nonrenormalizable. 
It is tempting to throw out the surface divergence as an artifact of idealized boundary condi- 
tions and retain the finite part as the true Casimir energy as is often done in the literature 
(see references in [3]). However, this is not a physically valid renormalization procedure. It 
has been shown that these surface divergences cannot be removed via renormalization of any 
physical parameters of the theory pQElE]. I n the zeta- function regularization technique, these 
surface divergences do not appear because in effect they are renormalized to zercQ. In special 
cases, like the parallel plate geometry with infinite plates, this is not an issue because in the 
limit as A — ► oo the Casimir force is finite. However, in any realistic situation where the plates 
are of finite size, the Casimir force diverges in the limit A — > oo. 

Unambiguous Casimir calculations can be carried out with idealized boundary conditions in 
the apparatus called the Casimir piston. A few years ago, Cavalcanti [5] showed for the case of 
a two-dimensional (2+1) massless scalar field confined to a rectangular region with Dirichlet 
boundary conditions that the Casimir piston can resolve the issue of nonrenormalizable surface 
divergences that appear in Casimir calculations. A Casimir piston contains an interior and an 
exterior region and Cavalcanti showed explicitly that the surface cut-off terms of the interior 
and exterior regions canceled. He also showed that the Casimir force on the piston is always 
negative regardless of the ratios of the two sides of the rectangular region. This is in contrast to 
calculations that can yield positive Casimir forces in rectangular geometries when the surface 
cut-off terms are thrown out and no exterior region is considered (see references in [1]). 

In this paper, we use a multidimensional cut-off technique [6] to obtain exact expressions for the 
cut-off dependent (A-dependent) part of the Casimir energy for a d-dimensional parallelepiped 
region. In the limit A — > oo, these yield nonrenormalizable hypersurface divergences and we 
show explicitly that they cancel out in the Casimir piston scenario for any dimension. We 
then derive exact expressions for the Casimir force on the piston in any dimension d and use 
the invariance of the vacuum energy under permutations of lengths to derive an alternative 

1 Like dimensional regularization, zeta-function regularization goes beyond pure regularization and does some 
renormalization. 
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expression. When the plate separation a is large, an otherwise long computation using the 
first expression becomes trivial using the alternative expression and vice versa when a is small 
(there is a useful a — > 1/a duality). As in two and three dimensions, the Casimir force on the 
piston is attractive (negative) in any spatial dimension d. 

For the three-dimensional Casimir piston with massless scalar fields obeying Dirichlet and 
Neumann boundary conditions approximate results were first obtained in [7J for small plate 
separation. Exact results for arbitrary plate separation were then obtained for the Dirichlet 
case in {£]. Exact results for the 3D Neumann (as well as Dirichlet) case were recently obtained 
via an optical path technique [12J. In [12j . arbitrary cross sections, temperature and free energy 
were also studied. We apply our d-dimensional formulas (both expressions) to the 2D and 3D 
Neumann cases. The first 3D expression looks similar in form to the one recently derived in 
|12j and is in numerical agreement with it. The alternative 3D expression converges quickly 
when a is large and though it is quite different in form compared to the first expression or 
the one found in [12] it is in numerical agreement with both of them. The 2D Neumann 
results are new and bring a completeness to the original work of Cavalcanti [5] where the 
2D Dirichlet case was considered. Before discussing the literature on Casimir pistons for the 
electromagnetic case, it is worth noting that the use of massless scalar fields in Casimir studies 
goes beyond theoretical interest and has direct application to physical systems such as Bose- 
Einstein condensates [20|. [2T1 [22] . Higher-dimensional scalar field Casimir calculations have 
also been carried out in the context of 6D supergravity theories [23J. 

For perfect-conductor conditions, the Casimir piston for the electromagnetic field in a three- 
dimensional rectangular cavity was studied in [7J and the Casimir force on the piston was 
found to be attractive in contrast to results without exterior region where the force could be 
positive. This was then generalized further in Refs. [9]- [12] where the temperature and free 
energy dependence was studied. It was shown in [13] that the Casimir force between two 
bodies related by reflection is always attractive, independent of the exact form of the bodies 
or dielectric properties and this was generalized further in [14]. It has also been shown that 
Casimir piston scenarios can yield repulsive forces. The Casimir piston for a weakly reflecting 
dielectric was considered in [15J and it was shown that though attraction occurred for small 
plate separation, this could switch to repulsion for sufficiently large separation. However, 
the force remained attractive for all plate separations if the material was thick enough, in 
agreement with the results in [7] . Two preprints [16} [T7] also discuss scenarios where repulsive 
Casimir forces in pistons can be achieved. Recently, two independent groups have developed 
techniques for calculating Casimir forces between arbitrary compact objects and have applied 
the results to the case of two spherical bodies at a distance [181 02] • 
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2 Cancellation of hypersurface divergences in the d-dimensional 
Casimir piston 



The expression for the vacuum energy E regularized using a multidimensional cut-off technique 
[6] divides naturally into two parts: a finite part Eq and a cut-off dependent part E{A) which 
diverges as the cut-off A tends to infinity. The expressions for E(A) and Eq are derived in 
appendix A and are written in a compact fashion with the help of the ordered symbol k 
which was introduced in [8| and is defined below after the expressions for E(A) are written. 
Our goal in this section is to show that for the Casimir piston scenario, the hypersurface 
divergences of the interior and exterior regions of the piston cancel out for any dimension d. 
By "cancel out" we do not mean that the cut-off dependent part of the Casimir energy is zero 
but that it is independent of the plate separation a so that the Casimir force on the piston has 
no cut-off dependence. In the next section we focus on the finite part Eq and obtain explicit 
expressions for the Casimir force on the piston in any dimension. We work in units where 
h = c=l. 

The regularized vacuum energy for massless scalar fields in a <i-dimensional box with sides of 
arbitrary lengths L\, L2, obeying Dirichlet (D) boundary conditions is given by (|A.25|) 
(as A — ► 00) 

E d = E 0d +E d (A) (2.1) 

where Eq d is the finite part for the Dirichlet case and Ed(A) is the cut-off dependent part 
given by (|A.13|) : 



1 d m 
m=l i=l 

1 d 

= ^TT E (~ l ) d+m m 2 ™ ^ V{m ^) Am+1 ZL,k m toi • • • LtJ. 

m=l 



(2.2) 



For the case of Neumann (N) boundary conditions the regularized vacuum energy is given by 
^A~2^ (as A -> 00) 

E n = Eq n +E n {A) (2.3) 
where Eq n is the finite part and Ej^(A) is the cut-off dependent part given by (|A.14[) : 

1 d m 

e n (A) es E m2-^T(m±l) A^ ^,..^11 ^ 

d m=1 i=1 (2.4) 

= ^rE m2m ^ r (^) Am+1 • • • L k J- 

m=l 
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In this section, the expressions for E$ D and Eq n are not needed. 

There is an implicit summation over the integers ki in (|2.2|) and (|2.4p . The ordered symbol 

id 



^,..,,,„ IB] is defined by 



d jl if A;i<fc 2 <...< A; m ; 1 < k m < d 

I U otherwise . 

The ordered symbol ensures that the implicit sum over the ki is over all distinct sets {k±, . . . , k m }, 
where the ki are integers that can run from 1 to d inclusively under the constraint that k\ < 
k 2 < • • • < k m . The superscript d specifies the maximum value of k m . For example, if m = 2 and 
d = 3 then £j£ fc m = fc 2 an d the non-zero terms are £1,2 , £ 1,3 and £ 2,3- This means the sum- 
mation is over {k\,k 2 } = (1, 2), (1,3) and (2, 3) so that fe ifc 2 = L\ L 2 + L\ L 3 + L 2 L 3 . 

2.1 Cancellation in three dimensions 

Before showing how the cut-off dependent hypersurface divergences in the d-dimensional Casimir 
piston cancel, we consider the case of three spatial dimensions first. Three dimensions allows 
us to make the first non-trivial use of the <i-dimensional cut-off expressions (|2.2p and (|2.4p 
and to illustrate in a transparent fashion how the cancellation occurs. This paves the way to 
follow the cancellation in (i-dimensions in the next subsection. The cut-off expressions in three 
dimensions and their cancellation in the piston scenario are in agreement with the work in 
[7J H2] and this provides an independent confirmation of our general formulas. 

In three dimensions, the Dirichlet and Neumann cut-off expressions Ejj(A) and E^(A), are 
obtained by substituting d = 3 in equations (|2.2p and (|2.4p : 



Ed(A) = ±J2 (-l) 3+m m2 m ^T(^)A m + 1 £f lr ., fem (L kl . . . L k J 

m=l (2.6) 

7T 7T 2 37T 2 

= - A 2 (Li + L 2 + L 3 ) - ^A 3 (Li L 2 + L X L 3 + L 2 L 3 ) + —A 4 L x L 2 L 3 
En {A) = ±fl e kl ,., km (L kl ...L k Jm2 m n^ T(^) A m+1 

Z m=l (2.7) 
7T 7T 2 37T 2 

= - A 2 (Li + L 2 + L 3 ) + y a3 ( L i L 2 + Li is + ^2 ^s) + —A 4 ii ^2 is • 

Except for a trivial redefinition of A — > A/7r, the above cut-off expressions in three dimensions 
are in agreement with those derived in |12j . The A 4 term appearing in (12. 6p and (|2.7p is 
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Figure 1: Casimir piston in three dimensions. 



multiplied by the volume L\ L2 L3 of the box and represent the volume divergence of the con- 
tinuum. This volume term poses no divergence problems since it must be subtracted to obtain 
the Casimir energy (defined as the difference between the vacuum energy with boundaries and 
the bulk vacuum energy of the continuum). In other words, it can be renormalized to zero via 
a constant counterterm in the Hamiltonian. The remaining A 2 and A 3 terms are proportional 
to the perimeter and surface area respectively (we refer to both as surface divergences for sim- 
plicity). In contrast to the volume divergence, there is no physical justification for subtracting 
out the surface divergences. In other words, they cannot be renormalized to zero. 

For the Casimir piston, the Casimir energy is obtained by adding the vacuum energy of the 
interior region I and exterior region II (see fig. 1). To obtain the cut-off dependence for the 
Casimir energy we therefore add the cut-off terms (the surface divergences) in regions I and 
II. Let the plate separation be a. In region I, the three lengths are L\ = a, Li = b and L3 = c 
whereas in region II the three lengths are L\ = s — a, L2 = b and L3 = c. Note that in region 
I, L\ comes with + a whereas in region II, L\ comes with the opposite sign — a. For Dirichlet 
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boundary conditions we obtain 



E Dpiston (A) = E Dl (A) + E D2 (A) 



_2 



= -A 2 (a + b + c) A 3 (ab + ac + be) 

8 4 

vr vr 2 ( 2 - 8 ) 

+ -A 2 (s -a + b + c) - ~^A 3 ((s - a) b + (s - a) c + be) 

_ 2 

= -A 2 (s + 2b + 2c) - ^-A 3 (sb + sc + 2bc) 
8 4 ; 



For Neumann boundary conditions we obtain 

E Nptston (A) = E Nl (A) + E N . 2 (A) 



= -A 2 (a + b + c) + — A 3 (ab + ac + bc) 
8 4 v ' 

+ -A 2 (s -a + b + c) + — A 3 ((s - a) b + (s - a) c + be) 

2 

= -A 2 (s + 2b + 2c) + ^-A 3 (sb + sc + 2bc) . 
8 4 



(2.9) 



Both Ejj piston (A) and E^ piston (A) have no dependence on the plate separation a. This is due 
to a cancellation that has occurred between region I and II. The Casimir force on the piston 
has therefore no dependence on the cut-off A (since the partial derivative with respect to a of 
E D piston {A) and E Npiston (A) is zero). 

2.2 Cancellation in d dimensions 

In a ci-dimensional Casimir piston, the piston has d—1 dimensions and divides again the volume 
into two regions: an interior region I and exterior region II. Without loss of generality, the 
direction in which the piston moves is chosen to be along the L\ direction so that region I and 
II share the same lengths except for L\. It is therefore convenient to write the Dirichlet and 
Neumann cut-off expressions (12. 2\ and (|2.4|) as a sum of two terms: one that includes L\ and 
another which is independent of L\ i.e. 

_. d mm 

e d{A) = ^ ^(-l)^m2«'^r(f) A-+ 1 ( Ll tU^H L ki +et,„^n L * 

m=l i=2 fc i^l i=l 



1 d mm 

__ /M ( Ll ^f, k2: .. :km ]jL ki +^ >km l[L ki 

m=l i=2 fc i^ 1 i=l 



(2.10) 
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where f(m) = m2 m it™^ Y(^^-) A m+1 . The Neumann cut-off expression is given by 

^ d m m 

^(A) - e /m ( l i ^....^n ^ +^,.^n • ( 2 - n ) 

m=l i=2 fci^ 1 i=l 

Let the plate separation be a. In region I, L± = a and in region II, L\ = s — a (the piston 
splits the length s into a and s — a along the L\ direction) (see fig. 1). To obtain the cut-off 
dependence for the <i-dimensional Casimir piston we need to add the contributions from regions 
I and II: 

E Dpiston (A-) = E Dl (A) + E D2 (A) 

^ d mm 

m=l i=2 fc i^l i=l 

. d mm ^ 2^2) 

+^E(- i )^/(H((«-a)^,* a ,., fcm n^+ft,.,^ii L * ( ) 

m=l i=2 fc i^l i=l 

^ d m m 

m=l i=2 ^i^ 1 i=l 

^iv piston (A) = ^ 1 (A) + £; i v 2 (A) 

^ d m m 

= ^ /(m) (° ^M,..,k m II Lfc * + £fci,.-,km II Lfe ») 

m=l i=2 ^i^ 1 i=l 

1 ^ m m (2 13) 

^PT E /M(( S - a ) Cl.fca,..,^!! Lfc * + ^fel,--,femll 

m=l i=2 i=l 

^ d m m 

m=l i=2 i=l 

The cut-off expressions for the piston, Eo viston {A) and E^ piston (K), have no dependence on the 
plate separation a. Their derivatives with respect to a is zero which implies that the Casimir 
force on the piston has no cut-off dependence in any dimension d. The hypersurface divergences 
have cancelled out in all dimensions in the Casimir piston scenario. 



+ 2 <m 



3 Casimir force formulas in the d-dimensional Casimir piston 

Having proved the cancellation of hypersurface divergences in the <i-dimensional Casimir piston, 
we now focus on the finite (A-independent) part of the Casimir energy. The finite part is 
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conveniently expressed as a sum of two terms: an analytical term composed of a finite sum 
over Riemann zeta and gamma functions and a remainder term Rj composed of infinite sums 
over modified Bessel functions (though convergence is reached after summing a few terms). 
In appendix A we derive exact expressions for the finite part Eq n and Eq d of the Casimir 
energy in d dimensions for Neumann and Dirichlet boundary conditions respectively. In this 
section, we state these expressions and use them to obtain the Neumann and Dirichlet Casimir 
force -F/v and Fjj for the <i-dimensional Casimir piston. In appendix B we develop alternative 
expressions for the Casimir force (see discussion at the end of this section). 

The finite part of the d-dimensional Casimir energy for Neumann and Dirichlet boundary 
conditions, Eq n and Eq d , is given by (|A.22j) and (|A.26l) respectively: 



_ d m—1 j j 

^tEE irw 1 ^^ & + 2 > + R *) ^ 



E, 



On 



m=l j=0 

where the remainder R^. is given by (|A,23D 



(L m y +1 



Rn, 



E E' 



2 n^r K i+1 



n=l 



-oo 



TT 



2tt n 



? ife) 2 + --- + (^x^) 2 ) 



+ 



4 



(3.15) 



and 



Eor 



TT 



d-1 



2 d+ 



1 Lk -'- Lk ^nm^aj + 2) + R Dl 



j=0 



where i?n is given by (|A.27P : 



(3.16) 



« DJ = E E 



2.f K l±l 



TT 



=— oo 
=1,-J 



2tt n 



H + 



(I 



1 L d 



+ ■ 



(3.17) 



The prime on the sum in (I3.15|) and (I3.17|) means that the case when all ^'s are simultaneously 
zero (l\ = £2 = ■ ■ ■ = lj = 0) is to be excluded. There is an implicit summation over the 
kiS via the ordered symbol ^k lt ..,kj defined in (12. 5p . R^. and ifo. do not depend only on j 
but are also a function of the ratios of lengths, for example i?jv, = RNjiL^/Lm, . . . , LkjL m ). 

For j = 0, 

and Lfc. are defined to be unity so that 



Therefore, the implicit summation over the k^s applies also to Rnj and RDj 



Rm 3 and ifo. are defined to be zero and ^ 
(L kl . . . L kj )/{L d y+ 1 = 1/L d for j = 0. 
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To obtain the ci-dimensional Casimir energy in the piston scenario we need to sum the contri- 
butions from region I and region II. In region I, let the length of the sides of the (i-dimensional 
parallelepiped region be a\, 0,2, ■ ■ ■ , a^-i, a where a is the plate separation. In region I, we label 
the lengths Lj such that L\ = a%, L2 = «2, etc. with = a {L& is equal to the plate separa- 
tion). In region II, the length of the sides are the same as in region I except for the length a 
which is replaced by the length s — a. The d lengths are s — a, a%, 0,2, ■ ■ ■ , a<i—l- For region II, 
we choose to label the lengths L, L such that L\ = s — a, L2 = ai, L3 = 02, . . . , = ad-i- To 
calculate the Casimir force, we only need to keep terms in the Casimir energy that depend on 
the plate separation a: in region I, this means keeping terms with = a and in region II this 
means keeping terms with L\ = s — a. 

In region I, the a-dependent Casimir energy for Neumann boundary conditions is obtained by 
setting m = d so that L m = = a and setting L^. = a^. in f)3. 14|) : 



E o Nl (*) = -^T ECU ^JW^ rf CO' + 2) + Ri Nj ) (3.18) 

3=0 



with Ri N . given by 



Rl Nj =1^1^ — 1+1 • (3 - 19) 



•j t-^i t < -ft 

n=l ti=— 00 
i=l,. 



1 ^ )2 + ... + (^.^ ) 



1 



A word on notation: the Roman numerals I and II will denote region I and II respectively 
while j will be denoted via Arabic numerals 1,2,3 etc. e.g. Rj N means the remainder(R) for 
Neumann(N) in region I with j = 1. 

The finite part of the Casimir energy for Dirichlet boundary conditions in region I is obtained 
by setting L& = a and L^. = in (|3.16p : 

e« d , = 4» E(-D" + ' f£U (r(¥) ^ CO + 2) + n lDi ) (3.20) 

where Rj D is given by ()3.17p : 



^,=E E — m • ( 3 - 21 ) 



7T 

n=l li=—oo 
i=l,... ,3 



(4^)2 + ... + (/ .fSL )a 



4 



In region II, only terms with L\ = s — a contribute to the a-dependent Casimir energy. We 
therefore consider only the cases when k\ is equal to 1 so that = s — a. The rest of the 
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lengths (J > 1) are given by L^. = cu-i so that = 01, L3 = a^-, ■ ■ ■ ,L m = a m -\. We are 
interested in an exterior of infinite length so that the Casimir force in region II is calculated in 
the limit s— >oo. For Neumann boundary conditions, the case (j = 0, m = 1) in (|3.14p can be 
omitted because L m = L\ = s — a appears in the denominator and yields a zero Casimir force 
in the limit of s — > 00. The cases j = 0, m > 1 can also be dropped because they do not yield 
any terms with L\ (with j = 0, the numerator is equal to unity and with m > 1, L m 7^ L\). 
In region II, the lower limit of the sums in (|3.14p therefore start at j = 1 and m = 2 yielding 
the following Neumann energy (a-dependent part) : 

E »»n («o = -^ee 2d ~ m e&u ( r (^) ^ co- + v + ^ ) 

(3.22) 

with 

00 00 i±l ( 2vr n . /(£ 1 ^2_)2 + . . . + (Z ^i)2 ) 

^ =z^z^ — m — • (3 - 23) 

For Dirichlet boundary conditions, the case j = can be dropped from (|3.16|) because it does 
not yield any terms with L\. The sum in (|3.16p therefore starts at j = 1 and we set = s — a 
and L fc = a k ,_ 1 , with L d = a d -i: 

J 3 

*« % w = D-d"^ «ti "^p" ( r <^> cw + 2) + r„ d ) (3.24) 

with 

=EET — — • (3 - 25) 

It is now straightforward to calculate the Casimir forces in each region. The Casimir force 
contribution from region I for Neumann is 

Fn ' ~dW = -^l^Zk u ..,k j ( a ki--- a kj)-]T2 r { V )t 2 CU + 2) ~— (3.26) 

3=0 

where Rj N is the remainder contribution given by 

vr d ~ 1 (afc a/e ) 
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The corresponding formula for Dirichlet are 

F "< = "ft"* = 2^ D-D*' « £■* («* ■ ■ ■ ^ r <^) + 2 > - 

(3.28) 

with given by 

The Casimir force from the exterior region II is obtained in the limit when s tends to infinity. 
For the Neumann case one obtains 

o d m-1 a „ 

rra=2j=l 

d m-1 a a 

'tee ^-«csu t":;,r «»„,<<. = « 



2^^^ (a m _iV' +1 

m=3j=2 v ; 



(3.30) 



where we used the resultH — lim s _^oo — — [(s — a) Rjx n .] = Rii n .{£i = 0). Rn N .(£i = 0) means 

(J (X 3 3 3 

Rji N , evaluated with i\ = 0. Note that the product a k _ x . . . a fc j that appears in (|3.3U|) is 
identically equal to one for j = l so that ^TV 1 u. {fl k _ 1 . . . a fc ._ 1 )/(a m _i)- 7+1 = l/(a m _i) 2 for 
3 = 1- 



For the Dirichlet case one obtains 



d-l 



F °- - ~ Os. = 5£rB-i>*" ig^T r(¥) co + n 



3=1 

V sA, ^d+n^d-1 %-l 



(3.31) 



2 -lim s _ 00 — [(s-a) -R/jjv 1 = lim^-nx, -lim^oo (s-a) — . The first term yields i?jj JV (li = 0) 

since the modified Bessel functions decrease to zero exponentially in the limit s — > oo except when = (as 
there is no s dependence when l\ — 0) . The second term is zero because the derivative of the modified Bessel 
functions with respect to a decrease exponentially to zero in the limit s— >oo when i\ 7^ 0. When £i = 0, Rn N . 
no longer has any dependence on a so that its derivative is zero identically. 
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The Casimir force Fry and Fu on the piston for Neumann and Dirichlet respectively is finally 
obtained by adding contributions from both region I and II: 

F N = F Nl +F Nn ; F D = F Dl + F Da . (3.32) 

Eq. (f332D together with (f3T26l) and (1330]) for F Nl and F Nn , and <^28\\ and (^BTT) for F Dl and 
Fpjj respectively constitute our final result for the Casimir force on the piston in d dimensions. 

The modified Bessel functions Kj+i that appear in Ri N or Ri D (Eqs. (|3.19j) and (|3.2ip ) 

2 3 3 

converge exponentially fast if the plate separation a is the smallest of the d lengths because 
the ratios a^./a that appear in the argument of the Bessel functions are then greater than or 
equal to 1. Only a few terms need to be summed to reach high accuracy and the result is also 
small in magnitude. This is why we can call R a remainder. However, Ri N and Ri D can 
converge slowly and be large if a is larger than the other lengths. In particular, the large a limit 
when a^/a« 1 would require a very large number of terms to be summed before convergence 
is reached. By making use of the invariance of the vacuum energy under permutation of 
lengths, we derive in appendix B alternative expressions that are more convenient to use when 
the plate separation a is large. For Neumann and Dirichlet they are given by (jB.lOh and (IB.13j) 
respectively: 

and 

F£ lt = -^-Rgf^O) (3.34) 

where R?J(ti^0) is given by ([R6]) with (|BTTT|> and Rfflii^O) is given by flRSj) with (lBl4|> . 
The above compact formulas are applied in the next section in two and three dimensions where 
one can see explicitly how they are used. 

The ratio of lengths that appear in the argument of the modified Bessel functions in (|B.11|) 
and (|B.14h have a in the numerator (a/a^) in contrast to our original expressions (with ratio 
aki/a). We have a useful a — * 1/a duality: when a is large a long computation with the original 
expressions can be trivial using the alternative expressions and vice versa when o is small. The 
invariance of the vacuum energy under permutations of the d lengths was used to derive the 
alternative expressions and the duality can be traced to this symmetry. Note that regardless 
of the size of the plate separation a, we would want to label the other d— 1 lengths such that 
o-i > 0,2 > as > . . . > a^_i to reach the quickest convergence. 

The Casimir force on the piston is negative (attractive) in all dimensions for both Neumann 
and Dirichlet boundary conditions and ranges from — oo (in the limit a^O) to (in the limit 
a — > oo ). The limit as a — > is easily determined using the original expressions (|3.32p . In 
the limit a -> 0, F Nl and F Dl given by (EHS]) and (l33H|) tend to -l/a d+1 ( dR lN /da and 
dRj D j da tend to zero) . F^ n and Fd u have no dependence on a so that Fry = F/v 7 + Fn xi 
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and Fd = Fjjj + -Fd 7J tend towards — 1/a +1 and hence — oo as a— >0. To determine the limit 
as a— > oo, it is easiest to use the alternative expressions. As already discussed at the end of 
appendix B, F^ lt and F£ lt given by (|3.33[) and (|3.34p tend to zero in that limit because the 
modified Bessel functions that appear in Rfjf^^O) (Eqs. ([B~l)|) and (lETlj) ) and R^(h^0) 
(Eqs. (|B.8p and ()B.14|) ) decrease exponentially fast to zero as a— >oo (since ^17^0, when a -^-oo 
the argument of the Bessel functions tend to infinity) . The rapid decrease to zero can be seen 
in the two and three-dimensional plots of the next section (fig. 2 and fig. 3). 

4 Application: two and three-dimensional Casimir piston for 
Neumann boundary conditions 

Exact results for massless scalar fields in the two and three-dimensional Casimir piston for 
Dirichlet boundary conditions were first obtained in Refs. [5j [8] and recently exact results for 
3D Neumann (as well as Dirichlet) were obtained in [12]. We apply our d-dimensional formulas 
(both expressions) to the two and three-dimensional Casimir piston with Neumann boundary 
conditions. The 2D Neumann results are new and fill a gap in the literature. For 3D Neumann, 
our first expression looks similar in form to the one recently derived in |12j and is in numerical 
agreement with it providing an independent confirmation of our results. Our alternative 3D 
Neumann expression looks quite different in form from the first expression and yields the same 
numerical results but is more useful (converges more quickly) at large plate separation a. 



4.1 Two dimensions 

In d dimensions the lengths of the parallelepiped are a, a±, 02, ■ ■ ■ , a^-i with a being the plate 
separation. In two dimensions the lengths are then a and a\ (we set a± = b so that the geometry 
is an a x b rectangle) . The Casimir force contribution from region I is obtained by setting d = 2 
in (f3T26l) : 

F »> = -I Eel,,,, («* ■••«*,) ^ r ( ^) ^ co' + 2^ dRl » 



ai+ 2 K 2 ' ^ ; da 

3=0 



8 a 2 '™' w ' 2 

n=l 1=1 



(4.35) 



= -- — 5--- L ^7 + -7^> > t -Ki(2n nib a) 
48 a 2 8vra 3 2 da ^ ^ i a v 1 ! 

n=l 1=1 

where Ri N is obtained from ()3.27|) : 

. 1 00 00 _. 

R lN = ~- 2 Ri Nl (b/a) = -^EEl- K ^ ni b M ■ ( 436 ) 
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Figure 2: Neumann Casimir force F versus a/b where a is the plate separation and b the 
second side of a two-dimensional rectangular region. The force is in units of 1/b 2 . The force is 
negative with magnitude decreasing quickly to zero as a/b increases. 



Ri N (b/a) is obtained from (|3. 19|) and means Ri N is a function of b/a. 

The Casimir force contribution from region II is obtained by setting d = 2 in (|3.3(J|) . In the 
first double sum, there is only the term j = l,m = 2 to consider. The second double sum is 
zero (it is nonzero only starting at d = 3). We obtain the simple expression 



FNn 



C(3) 
16vr6 2 



(4.37) 



By summing F^j and F^ n we obtain the Casimir force Fn on the piston: 



F N 



7T 



((3)6 ((3) , 1 8 



48 a 2 8?ra 3 16vr6 2 2 5a 



(4.38) 



71=1 £=1 



As an example, the above expression yields F N = -0.1342935575/6 2 for the case of a square 
(a = b). 

An alternative expression Ffi lt for the Casimir force can be obtained via (I3.33p . For d = 2 we 
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obtain 



, oo oo (4.39) 



24a 2 26 ^ ^ t 

n=l 1=1 

where the prime on the modified Bessel function means partial derivative with respect to a i.e. 
K'(x) = jfoK(x). RfJ(£i^0) is obtained from the j = l, m = 2 term in (|R6l) : 

oo oo 
n=l i=l 

(^1^0) i s obtained from the j = l, m = 2 term in (IB.llh . 

The two expressions, (|4.38|) and (|4.39|) . yield the same value for the Casimir force on the piston 
and are valid for any values of a and b. However, computationally, expression (|4.38j) is better 
to use when a is small (i.e. a/b< 1), whereas (14.39)) is better to use when a is large (b/a< 1). 
This is the simplest case of the a — > 1/a duality that was discussed last section. 

The Casimir force on the piston is negative (attractive) and ranges from — oo (in the limit 
a— s-0) to (in the limit a^oo). A plot of F versus a/b (in units of 1/b 2 ) is shown in fig. 2. 



4.2 Three dimensions 

In d-dimensions we have the d lengths a, 01,02, • • • , dd-i with a the plate separation. In three 
dimensions the three lengths are then a, a\ and a<i- For the three-dimensional Casimir piston 
it has become customary to use a, b and c for the lengths and we therefore set 02 = b and 
a\ = c. The Casimir force contribution from region I is obtained by setting d = 3 in (13.261) 

2 

i=o ' (4.41) 

= _ C(3) (6 + c) 7T 

480 o 4 16?ra 3 96 a 2 /jv 
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where Rj N is given by (|3.27[) and f)3. 19j) : 



7T 



Ri N - 



c b 2 c & c 

3 ^ (c/a) + -2 ^ (6/o) + (c/6) + -3 Ri N2 (c/o, 6/c 



a 

00 00 



-yy" ? riKi(2vrn^c/a) + - KAI-k nib/ a) + % K^irnic/b) 
4 ^-^ £ I a a v 6 y 



n=l l=\ 



be 



n=l £i,£ 2 =— co 



(4.42) 



a J 



±1 
a 



+ 



•6\2 



-1 3/4 



i?/^ (c/a) means i?/^ is a function of c/a and the prime above the sum means that the case 
^1 = ^2 = is to be excluded from the sum. 

The Casimir force contribution from region II is obtained by setting d = 3 in (|3.30p : 

3 m— 1 



k2 ~ 1 "' fej - 1 r(^) 7r ^ c(j + 2) 



m=2 j=l 

n=l fcl 



where (13.23h with j = 2 and m = 3 was used to obtain 

( 

n=l ^=1 



(4.43) 



(4.44) 



The Casimir force Fn on the piston for Neumann boundary conditions in three dimensions is 
obtained by summing Fn x and Fn h ■ 



N 



ir 2 be C(3)(6 + c) vr 



480 a 4 16vra 3 

2 



96 a 2 R ' In 



7T 



00 00 



14406 3 167rVc 2 26 2 / 4b 3 ^ ^\£cJ 

n=l £=1 



(4.45) 



where Ri N is given by (|4.42p and R' In = dRi N /da . Note that the third term in square brackets 
in (|4,42p depends only on b and c and can be dropped when evaluating R^ . For the case of a 
cube (a = b = c), Eq. (|4.45[) yields Fjy = — 0.1380999/c 2 . As in two dimensions, the force Fn 
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Figure 3: 3D plot of Neumann Casimir force Fn versus a/c and b/c. The force is in units 
of 1/c 2 . The force is large and negative at small values of a/c and remains negative with its 
magnitude decreasing quickly to zero as a/c increases. The value of b/c shifts the magnitude 
of the force towards larger values as it increases. 



is negative and ranges from — oo (in the limit a^O) to (in the limit a^oo). A 3D plot of 
-F/V versus a/c and b/c (in units of 1/c 2 ) is shown in fig. 3. 

Our exact expression (|4.45j) looks similar in form to the one derived in [12] and is in numerical 
agreement with it. Moreover, in the small a limit, Rj is exponentially suppressed (exactly zero 
in the limit a— >0) and when b = c, the second row in (|4.45p yields 0.0429965/c 2 in agreement 
with the Neumann results in both [7] and |12j . This provides an independent confirmation of 
our results. 

An alternative expression F^ lt for the Casimir force can be readily obtained by substituting 
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d = 3 in (HEMP and using ([R6]) and (IBTTI) : 



n ait 

Af — — 



24 a 2 



7T 



9a 



9 a c 
~da~ IF 



24 a 2 



7T 



n=l£ 1= l£ 2 =-oo (£ia/6) 2 + (£ 2 c/6) 2 




(4.46) 



The prime above K denotes partial derivative with respect to a. Note that the sum over 
£2 includes £2 = (since the sum over £2 contains no prime above it). Equation (14.460 is 
our alternative expression for the exact Casimir force on the piston in three dimensions for 
Neumann boundary conditions. It is valid for any values of the lengths a, b and c and yields the 
same Casimir force as the original expression (|4.45p . However, in the argument of the modified 
Bessel functions the plate separation a now appears in the numerator making the alternative 
expression converge quickly for large a. A long computation with the original expression when 
a is large can converge exponentially fast with the alternative expression and vice versa when a 
is small. This is a nontrivial case of the a — ► 1/a duality already encountered in two dimensions 
and discussed in general in the last section. Note that we are free to label the base such that 
c > b to obtain the best convergence. 

5 Summary and discussion 

By applying a multidimensional cut-off technique we obtained expressions for the cut-off de- 
pendent part of the Casimir energy for parallelepiped geometries in any spatial dimension d 
and showed explicitly that nonrenormalizable hypersurface divergences cancel in the Casimir 
piston scenario in all dimensions. We then obtained exact expressions for the d-dimensional 
Casimir force on a piston for the case of massless scalar fields obeying Dirichlet and Neumann 
boundary conditions. As an example, we applied the d-dimensional formulas to the 2D and 
3D piston with Neumann boundary conditions. The two main features of the Casimir piston 
originally mentioned by Cavalcanti [5] for a 2D rectangular geometry, namely the cancellation 
of the surface divergences and the negative Casimir force on the piston, were shown to hold 
true in all dimensions d. We obtained two different expressions for the d-dimensional Casimir 
force. The Casimir energy is clearly invariant under permutations of the d lengths of the paral- 
lelepiped. This symmetry is trivial but its application is very useful: one can derive alternative 
expressions for the Casimir force that converge quickly compared to the original expressions 
when the plate separation a is large. 
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It would be interesting to generalize our ci-dimensional results to include arbitrary cross sections 
and thermal corrections. For scalar fields in three dimensions this has been recently considered 
in [12] via the optical path technique. The scalar field results were then used to obtain the 
electromagnetic (EM) Casimir energies with perfect metallic boundary conditions |12| . It would 
be worthwhile to see how the 3D alternative expressions for massless scalar fields derived here 
for Neumann and for Dirichlet elsewhere [8] can be modified to include thermal corrections. 
These could then be used to obtain 3D alternative expressions for the thermal corrections to 
the EM case. 



A Cut-off dependent and finite parts of the regularized vac- 
uum energy: periodic, Dirichlet and Neumann boundary 
conditions 

We consider a massless scalar field confined to a <i-dimensional parallelepiped region with ar- 
bitrary lengths L%, obeying periodic, Neumann and Dirichlet boundary conditions. Our 
goal is to include the cut-off dependent and finite parts of the vacuum energy regularized via 
a multidimensional cut-off technique [6J. This appendix naturally divides into two parts. We 
first consider periodic boundary conditions and make use of formulas found in section 2 and ap- 
pendix B of [6]. In particular, we determine explicitly the d-dimensional cut-off dependence in 
the expression for the regularized vacuum energy. In contrast to dimensional or zeta-function 
regularization, the multidimensional cut-off technique performs no renormalization. The sec- 
ond part consists of finding the regularized vacuum energy for the Dirichlet and Neumann 
cases. This is obtained by summing over the vacuum energy of the periodic case. 

The vacuum energy for periodic boundary conditions regularized using a cut-off A is given by 

m 



-.2 /„2 „2 



\ / 1 i i __d OO \ / 1 i I '"d 



rii=— oo rii=—oo 
i=l,...,d i=l,...,d 



oo ln\ oo oo . h\ n\ oo oo . h\ 

n\=—oo ri2=— ooni=— oo n^=— oo rii=— oo 

i=l,...,d-X 



J i 



d-l 
3=0 



(A.l) 
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x oo , / n 2 n\ rij 



where 

A,(X )S E= V^ + If+ +L >- (A.2) 

n=— oo rii =— oo 
i=l,...,j 

The prime on the sum over n means that n = is excluded from the sum. The function 
d\ Aj(A) can be expressed in the following form [B] (in the limit A — > 0) 



d x Aj(X) = ^^L( 2^ d x ^2 r e- X '^ + - +x hx 1 ...dx j + R j ) (A.3) 

\ L j+l) J v ri=l Jo 

where A' = \/Lj + \ and L\. . . Lj is a product of lengths i.e. Ylj =1 Li. This product is defined 
to be unity for the special case of j = 0. Rj is given by [6] 

oo oo J~H 

*i = £ £' 2( 2 jTr^2±i(ftrr^ Ll ) 2 + --- + ^^ 2 )- 

n=lt l =-ooTT[{l l L 1 ) 2 + --- + {l j L j ) 2 }— 2 
i=l,...,j 

(A.4) 

i?j starts at j = 1 (it is zero for j = 0). The functions Krj + iy 2 are modified Bessel functions 
and the prime on the sum means that the case where all the £iS are zero is excluded. Via the 
Euler-Maclaurin formula, the integral term in the round brackets in (|A.3|) can be decomposed 
into a cut-off dependent term (which diverges as A — > 0) and a finite term which is independent 
of the cut-off (see section 2 of [6]): 



, dxj 



= r ( ^)rfc(j + 2) + ^2^r ( ^)-^2^Mr ( ^) 

= r(^)rf C (j + 2) + ^2*^ r(^) (v)^ 1 - ^^ir* rffi (v) j+2 . 

(A.5) 

Substituting (|A.5[) into (|A.3[) . the regularized vacuum energy (jA.ip for periodic boundary 
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conditions is given by ( as A — > 0) 



d-l r r 



£5 



r(2+2 )7r ^ c(j + 2) + jRj 



+ 2 J ^ r(i±l) - ^ + 11 2^ vr§ r(^) (L j+1 )^ 



- E (^5^ ( r (^) ^ CO' + 2) + Rj) +(la... L d )j*z 2* r(^) 



^-^^-^^+...-i ?1 ^_ i?2 !^ + ... + ^ 2 ^^r(^)(L 1 ...i, / , 

6Z/i 2vr L| 90 Ljj L% L\ A d+1 V 2 ; 

(A.6) 



The notation £"„ (d, A) is a compact way of stating that the regularized vacuum energy 

f L 1 ,...,L d 

E p is a function of the dimension d, the cut-off parameter A and the lengths L\, Lj. 

The regularized vacuum energy for Dirichlet and Neumann boundary conditions can be ex- 
pressed as a sum over the periodic energies E p [8j [23] : 

E = E (± I)""" « fc. 4^.^ ('». ^ (A.7) 

m=l 

where the plus and negative signs correspond to the Neumann (N) and Dirichlet (D) cases 
respectively. £ fcm is called the ordered symbol and is defined by [2D] 

d jl if A;i<fc 2 <...< k m ; 1 < fc m < d 

I U otherwise . 

The ki's are positive integers that can run from 1 to a maximum value of d. The ordered symbol 
£fei fc m ensures that the implicit summation over the fcj's is over all distinct sets {fei, .., A; TO } 
under the constraint that k\ < k<i < ■ ■ ■ < k m . The superscript d specifies the maximum value of 
k m . E p (A, m) is obtained from (IA.6I) by replacing d by m and Lx by , L 2 by , 

etc. When substituted into (|A.7|) one obtains 

" 7 1 J= ° J+ ' (A.9) 

+ ^rE(± ^ • • • ^J^i 2m ^ r (^) 

m,=l 
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where Rj is now the function (|A.4|) with L\ replaced by L kl , L2 by Lk 2 , etc. It is convenient 
to define the function 



U fcl ,.., fc , +1 - ( r (^) ^ CO' + 2) + a,) . (A.10) 

Using (lA.lOj) and rearranging the limits on m and j, we can express (IA.9j) in the following 
compact form (we write out separately the Dirichlet and Neumann cases) 

Ed = ^EE ^) d+m ZL, k J 3kl ,.., k3+1 +Eo{K) (A.ll) 

j=0 m=j+l 

E * = ^L.,Uj kl ,.., kj+1 +E N (A) (A.12) 

J=0 m=j+l 

where the functions Ed (A) and Epj(A) are the cut-off dependent terms for the Dirichlet and 
Neumann cases respectively obtained from the second row in (|A.9|) (we now work with the 
cut-off A = 1/A instead of A so that the divergent limit A — ► is replaced by A — > 00 which is 
the more customary notation): 

1 d 

Ed(A) = ^Y.(- 1 ) d+m ZL.,k m (Lk l -..L k Jm2 m n^T(^)A m+1 (A.13) 

m=l 

E N (A) ee ^^L^fe-^J m2m ^ r ( ! f)A m+1 (A.14) 

m=l 

Note that in (|A.12p the limits on m and j in the double sum have been rearranged compared 
to (|A.9p . We can decompose k into a sum of two terms: k _ <2 + fc m - ^ n * ne 
first term, fc TO is set to its maximum value of d and the implicit sum is over the remaining fcj's 
with the maximum value of k m _i equal to d— 1 (hence the superscript d— 1). The second term 
contains the remaining implicit summation with the maximum value of k m equal to d— 1 (hence 
there is a superscript d— 1 in the second term as well). For the case m = d, the decomposition 
yields only one term £^ fcd = ^f" 1 fcd d + since can only be equal to d. 

With this decomposition the sum over m becomes 

E (± !)" +m = E (± !)" +m [C^-^ + ] • (A.15) 

m=j-\-l m=j+l 

There are two separate cases to evaluate above: the plus sign (the Neumann case) and the 
minus sign (the Dirichlet case). The Dirichlet case has already been calculated in [20J and the 
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result is 

m=j'+l 

For the Neumann case one obtains 
d 

£fci,..,fcm ~~ ^ki,..,kj,d ( ^fci,..,fej + i ^fei,..,fcj + i,d ) , , , 

m=j+l (A.17) 



Each pair of round brackets contains the sum of two terms which are equal. For example, 
consider the first pair of round brackets (^f" 1 k + Cf" 1 k+1 d) ■ ^he ^ ac * ^hat %+ 2 * s e< l ua l 
to (i in the second term is irrelevant since the summation over fj in (|A.12p ends at 

A; J+1 for a given j. Therefore, ^f" 1 fc d is equal to CJ?" 1 k . +1 - The same logic applies to the 
other pairs of round brackets. Equation (|A.17p reduces to a recursion relation 

d d-1 



m=j+l m=j+l 



d-2 

£fei,.!,fcj,d + 2 (^fi,.%j,d-l + 2 E • 

m=jr'+l 



(A.18) 



Applying the above recursion repeatedly (another d—j— 2 times) yields 

E E ^ ( a -!9) 

m=j+l m=j+l 

With (|A.19p . the double sum in (jA.12|) can be expressed as 

d— 1 d _ d to— 1 

— 7r \ ' \ nd—mtm—\ f ^ \ \ nd—mcm—l f /a 9fi\ 

j=0 m=j+l ** to=1 j=0 

The function fj k fc is given by (|A.10p with kj+i equal to m. Substituting (|A.20p into 
(|A.12p yields our final expression for the Neumann regularized vacuum energy 

E n =E 0n +E n (A) (A.21) 

where the finite part Eq n is given by 

d m— 1 

e 0n = ^rEE 2 "" m /// ^ . i ; (r(¥) ^ CO' + 2) + • (A.22) 

to=i i=o v m; 
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The function Rjy is given by (|A.4|) with L\ — > L^, Lj + \ — ► = L m : 

^ =Et^ 4^ • ( A - 2 3) 

n=1 fc~ [Wfe) 2 + --- + (^fe) 2 ] 4 

The Dirichlet case is obtained by substituting (|A.16P in (jA.lip yielding 

j'=o 

The function /j fc fc d is obtained by setting kj+i equal to d in (|A.10j) yielding the Dirichlet 
regularized vacuum energy 

E d = E 0d +E d {A) (A.25) 

where the finite part Eq d is given by 

^ = £t D-u** «£U w( r( ^» C(J + 2> + ■ (A - 26) 

The function Rdj is given by (|A.4j) with Li — > L^, Lj+\ — > Lfe- +1 = E^: 

R A W 2n^ ^(2^^^ + ... + ^^ ) 

Kd j = 2^ in ■ \ A - Z1 > 

For j = 0, Rn 3 and i?^ are defined to be zero and £ki,..,kj an d are defined to be unity. 

The final expressions for the regularized vacuum energy are (|A.21|) for the Neumann case and 
(|A.25P for the Dirichlet case with the cut-off dependent parts -Ed (A) and En (A) given by 
(|A.13P and (|A.14|) and the finite parts E 0n and E 0o given by (jA.22j) and (|A.26j) respectively. 

B Alternative expressions for the (/-dimensional Casimir piston 

In section 3 we derived expressions for the Casimir force on the piston. In this appendix 
we derive alternative expressions. This is accomplished by labeling the lengths Li in region 
I differently compared to section 3. The Casimir energy is invariant under permutation of 
lengths so a different labeling does not alter the value of the Casimir energy. However, the 
different labeling leads to an expression with a different form. 
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In section 3 we labeled the lengths Lj in region I as L\ = a%, L2 = 0,2, ■ ■ ■ , L^-i = a^-i with 
equal to the plate separation a. We now label L{ such that L\ = a, L2 = ai,L^ = 02, . . . , La = 
ad-i so that L\ is equal to the plate separation. We do not change the labeling in region II 
(i.e. L\ = s — a, L2 = a%, L3 = 02, etc.) so that we only need to obtain new formulas for region 
I. 

The expressions for the finite part of the Casimir energy are given by (|3.14p and (|3.16p for 
Neumann and Dirichlet respectively. We only keep the a-dependent terms and for region I this 
means keeping only those terms with L\ = a. For Neumann, (|3.14p divides into two sums: the 
j = 0, m = 1 term where L m = L\ = a appears in the denominator and all other terms where 
L\ appears in the numerator (this occurs when k\ = 1 so that = a for j > with the other 
lengths given by Lj, = cik -i)- This yields (with "alt" as superscript for "alternative") 



d m—1 



2 d 



1+1 Z^ Z^ z 5fc i: 



m=l j=0 



_i Lja ■ ■ ■ 

-*i (L m )i+i 



Y{l^)^C(3 + 2) + R N] 



7T 



d m—1 



24 a 2 d 



^ \ \ nd—m £m—l 

i+i Z^ Z^ ^iM,- 



a a. , . . . a, . 



m=2j=l 



kj (a m -i) j 



+1 



(B.1) 



with 



00 00 

E £' 

n=l £ i =—oo 
i=l,...,j 



, 2n^ K^nJih^ + .-. + ^^f 



TT 



j+l 



(B.2) 



For Dirichlet given by (|3,16p . the case j = does not yield any L\ terms so that it can be 
dropped. L\ = a appears only in the numerator via = a when k± = 1 (with the other 
lengths given by = a&._i). We obtain 



d-l 



2 d 



(B.3) 



with 



n=l ii=— 00 
i=l,...,j 



_^_)2 + ... + (*. ^i" 1 



r) 2 ) 



7T 



^) 2 + -" + 



.7 + 1 

4 



(B.4) 
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The alternative expression for the Casimir force in region I for Neumann is 



- + ^rE E *~<ZU rff) CO + 2) - ^£ 



24a 2 2 d +!^^ S1 ^--*i (a m __iV' +1 v 2 7 w ' 5a 

m=2j=l 



(B.5) 



where 



d to— 1 aa a 

D aZ< ^ r,d-m t-m-l k 2~ 1 k j~ 1 r>alt a\ 

R i n -- ¥ TT2^2^ 2 ZiM,~,ki < am % • ( B - 6 ) 

m=2j=l \ m ij 



The corresponding alternative expression for Dirichlet in region I is 



where 



d-i n n pmcdt (B.7) 

^ V( l) d +i fi^ 1 k2 ~ 1 '- k >- l T(i +2 \^ C(i I 2) Jc 
2W^ ( ij ^-^ (a^xp+i i(2j7F CU + 2j 9a 

d-1 

To obtain the Casimir force on the piston we need to add the contribution from region II: F^ n 
and Fjj n given by (|3.30p and (|3.3ip respectively. For Neumann we obtain 

, alt jp alt , T? n OU In 71 S^S^nd-m,m-l %-l ' " %-l 



jp alt _ r? alt , jp _ H ~_±K \ " \ " n d-m c m-l fc 2-! fcj- 1 p /« n N 

(B-9) 



where Rn N (£i = 0) is (13.23H evaluated at t\ = 0. The above can be reduced to a more compact 

expression by noticing that the l\ = contribution to —dRf l ^/da cancels out with the last 
term in (|B.9|) . The alternative expression for the Casimir force on the piston for Neumann 
boundary conditions reduces to 

< = -^-^<(^0)- (B-10) 

To evaluate dBTOl) we exclude t\ = in (|B.6|) so that Rf given by (|B.2p is evaluated without 
including l\ = i.e. 



OO CO oo 



4^ ^ 2 ± i (2vrn^ 1 ^ T ) 2 + ... + (^2-i) 2 



^(w=ee E — ;:r • (b.id 



n=l ^i=l £i=—oo 
i=2,...,j 



v 1 a m _i / V J a m _i / 
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Compared to (|B.2|) . there is no longer a prime on the sum over l{ and % starts with the integer 
2 instead of 1. For Dirichlet one obtains 

FB lt = < + F DlI = -- + ^ E(-l) d+ ^S,.., fcj ^'"y+i" 1 = 0) (B.12) 

where Ru D {l\ = 0) is (|3,25|) evaluated at i\ = 0. The above can also be reduced to a more 

compact expression since the £\ = contribution of —dRf l ^/da cancels out with the last term 
in (|B.12j) . The alternative expression for the Casimir force on the piston for Dirichlet boundary 
conditions reduces to the simple expression 

0R a ^ 

F alt = _fe ( ^ 0) (R13) 

oa 

To evaluate (|B,13P we exclude i\ = in (|B.8P so that Rfj^ given by (|B.4p is evaluated without 
4 = i.e. 

<;iw«) = EE E — ^^-^ ^r^- cb.14) 

""'-'fa;- [«i3fc) 2 + --- + ^S?) 2 ] 4 

Our alternative expression for the Casimir force on the piston for the Neumann case is F{f l 
which is given by (fB~T0|) together with (fB~6|) and (|BTTT]) . For the Dirichlet case the alternative 
expression is F£ lt which is given by (|B.13|) together with (]B.8[) and (|B. 14[) . It is now trivial to 
see that in the limit as the plate separation a tends to infinity that the Casimir force on the 
piston is zero since the modified Bessel functions and their derivatives that appear in F£ u and 
Fp lt via (IB.llj) and (1B.14|) decrease exponentially fast to zero as a tends to infinity. 
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